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Abstract. - We propose an interpretation of previous experiments and numerical experiments 
showing that, for a large class of systems, distributions of global quantities are similar to a 
distribution originally obtained for the magnetization in the 2D-XY model [1]. This approach, 
developed for the Ising model, is based on previous numerical observations [8]. We obtain an 
effective action using a perturbative method, which successfully describes the order parameter 
fluctuations near the phase transition. This leads to a direct link between the D-dimensional 
Ising model and the XY model in the same dimension, which appears to be a generic feature 
of many equilibrium critical systems and which is at the heart of the above observations. 
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Following a first observation by Bramwcll, Holdsworth and Pinton [1] , many studies report 
that the probability density functions (PDF), for spatially or temporarily averaged quantities 
in correlated equilibrium [2] and out-of-equilibrium [3-5] systems have a generic asymmetric 



form similar to the so-called BHP distribution. This distribution (see figure 1(b) I is obtained 
from the magnetic fluctuations in the two-dimensional (2D) XY model of magnetism, in the 
spin wave approximation, in the zero temperature limit [6]. This " superuniversality" is clearly 
incompatible with the notion of universality classes in critical phenomena and its ubiquitous 
presence in out-of-cquilibrium phenomena appears mysterious to say the least. In fact it is 
easy to find critical systems where the distribution is radically different and in most situations, 
some deviation from BHP distribution is apparent. Even in the case of the 2D-XY model itself 
it has been recently established that small temperature-dependent corrections exist [7,9,10]. 
There are therefore strict physical criteria associated with both the observation of this generic 
behaviour and with the deviations from it, just as in the case of Gaussian fluctuations through 
the application of the central limit theorem. In this paper we expose these criteria through 
the study of a well known and well controlled equilibrium system, the 2D Ising model, the 
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aim being to show microscopically how a XY-likc behaviour appears in the Ising model. 
Using numerical simulations for the 2D Ising model [2, 8] we established that there exists a 
range of temperatures T*(L), or applied field H*(L) close to the critical point( 1 ) (T C (L), H = 
0), where L is the system size, where the PDF is similar to the BHP function. At this 
specific temperature or field the magnetization shows intermittent behaviour, where coherent 
structures appear on intermediate time scales, in analogy with injected power fluctuations in 
enclosed turbulent flow [3,8]. This should be contrasted with the behaviour at the critical 
point, where the amplitude of the structures and the ensuing intermittency are cut off by the 
boundaries of the available phase space. The cross over from approximate supcruniversality 
to universality class dependence is related to this change of regime. We approach the critical 
point from the ordered phase, making a perturbation expansion about the ordered state. The 
calculation shows that there is indeed a quantitative similarity between the fluctuations for 
2D XY and Ising models up to this threshold and this is the origin of generic behaviour in 
the Ising system. Through this result we are able to make some precise statements about the 
criteria leading to such generic behaviour in more disparate and less well controlled systems. 
We consider the 2D classical Ising model on a square lattice of size L x L described, after the 
Hubbard-Stratronovitch transformation, by N = L 2 continuous variables <fii, leading to the 
partition function Z [11]: 



/detK\ 1/2 
Z cx \ j I V<p exp 
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(1) 



where K is the coupling matrix whose elements are Kij = K(2\8 rur . + ^2s^ r i< r i+ s )> an< ^ 
(3 = K = 1/T. The vectors S are the lattice unit vectors ±£ and ±y, and A is an arbitrary 
parameter introduced by the Hubbard-Stratonovitch transformation. It is chosen so that all 
the eigenvalues of K are positive, and the mean-field critical temperature T c mf = 2(D + A) is 
usually defined with A = 0. The local spin cr, is thus mapped onto the local magnetization 

mi = tanh [J2f=i Kij<f>j + P H ) . 



Generalized Ginsburg criterion for fluctuations at large length scale. - The usual criterion 
used in critical phenomena to discuss the validity of a perturbative approach to the fluctuations 
is the Ginsburg criterion. This is defined such that the ratio of the two-point correlation 
function, averaged up to the correlation length £ and the square of the magnetization, averaged 
on the same scale, be small. For the 2D Ising model it is well known that this ratio is not 
small compare to unity, and therefore a perturbative approach can not capture the physical 
behaviour of the model up to this length scale. In our particular case however we are interested 
in the fluctuations at a large scale I, not equal to £. It is hence natural to define a generalized 
Ginsburg criterion by the ratio (for D< 4): 

Ri= j° :L, ■ (2) 



f Q (m(r)) 2 dr 



Here, G(r) is the two-point correlation function which is roughly equal to l/r D ~~ 2+v for r < £ 
and for r > £. (m(r)) is the local magnetization, which scales like ' . The second and 
approximate equality is obtained using the scaling hypothesis and is true independently of 

( )T* was defined in ref. [8] as the temperature where the kurtosis of the distribution most closely approx- 
imates to that for the BHP function. As usual for finite size systems, we define the critical temperature the 
one corresponding to the maximum value of the susceptibility. 
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the value of the critical exponents. The traditional Ginsburg criterion is recovered for I ~ £ 
and it is unsatisfied for all dimension D< 4 (A ~ 1). The main result of ref. [8] is that along 
the locus of temperatures T*(L), while the correlation length diverges with L as expected 
for a critical phenomenon, its amplitude is small compared to L: f/L ~ 0.03. Therefore we 
have Rl{T*(L)) <C 1 and we expect that the fluctuations at the integral scale I = L could be 
described, to an excellent approximation, by a pcrturbative development. Such an approach 
cannot of course give correct critical exponents but this is consistent with the observation that 
such behaviour is, to a very good approximation, independent of the exponents at hand [2,9]. 
From the point of view of renormalisation, we are expanding about the zero temperature 
fixed point, which is Gaussian. The major contribution to the fluctuations will come from 
length scales of the order of the correlation length. From the microscopic scale up to £ the 
effective action stays close to the non-trivial fixed point, only to cross over towards the zero 
temperature fixed point for I ~ L. The calculation approximates the behaviour at this scale 
with contributions estimated from the zero temperature sink. 

To make the pcrturbative development we define the "fast" , fluctuating variables 6i = 4>i — 0o 
which we separate from the "slow" variable </>o = Si <t>i/N. Here the terms "fast" and "slow" 
are taken from the dynamical point of view as we will see in detail below. We take the 9i to 
be small. The effective action, expanded to lowest order in the 0i and Fourier transformed 
can be separated into two parts: 

Z cx / d</> [ d0,d0_ g exp -S[<h,0 q ], (3) 
S[<h,0 q ] = S (c/> ) + E <%% -v G q 1 = 1 - ^-T 2 )K q , 

5o(0 o ) = N [K(D + \)<f>l - Iogcosh(2if (£> + A)0 O + /3H)] , 

where Si is the minimal set of Fourier modes q = (2Tvp/L,2ivq/L) ^ (0,0) that, in addition 
to the operation q = (2n(L — p)/ L,2tt(L — q)/L), fills up the entire Brillouin zone except 
the zero mode. This prevents us to count twice the quantities 9± q . We also define T = 
tanh[2if (D + A)0o + (3H] and K q = 2K(X + cosq x + cosq y ). The total magnetization is then 
given by m tot ~ T - 2T(1 - T 2 ) Z qeSl K q q 0- q /N. 

Dynamical approach. - To this lowest order, the action © looks similar to that of the 
2D-XY model with a propagator that is a function of H and (f>Q. Setting 0o constant would 
make it truly XYdike, with a massive propagator. However, at zero field this is not the case as 
we have here a finite size system. Hence, in dealing with the effective action we have to take 
into account the fact that there is no rigorous symmetry breaking and that the equilibrium, low 
temperature magnetization is strictly zero, in zero field. For this reason we choose a dynamical 
approach separating time scales for fluctuations about a local free energy minimum from those 
for a passage from one local minimum to the other. This separation of scales defines the fast 
modes and slow modes of evolution. At low temperature, or in finite field, the slow, or ergodic 
time scale will be outside the numerical or experimental observation time scale and we expect 
fluctuations around a single minimum. As the critical point is approached, we expect this 
separation of scales to be no longer possible with the result that the symmetry is restored. 
We can now follow the Langevin dynamics from the action J2J. Defining the fields <jy q (resp. 
(p^) as Re (9 q ) (resp. Im (0 q )), we obtain the following Langevin equations for the fast and 
slow degrees of freedom: 

#>(t) = -2G-^(t) +„£»>(*), Mt) = - + V o(t), (4) 
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(a) Comparison at T < T c and H = 



(b) Comparison at T = T* and H = 



Fig. 1 - Numerical test of the Langevin dynamics, for A = 2.5 and L = 65. (a) The data (x) are 
Monte Carlo simulations at T = 1.7 and L = 64. The data (•) are obtained by Langevin dynamics from 
equation (E) at T dyn = T c mf x 0.74. The plain curve corresponds to the distribution of the Gaussian 
model (L = 64) with mass M = 3.5 x 1CT 1 [14]. (b) The plain curve is the BHP distribution [6]. 
The data (x) are Monte Carlo simulations for L = 64 system at T*(L) = T c x 0.94 [8]. The data 
(•) are from the Langevin dynamics JIJ, at T*(L, A) = T c mf x 0.95. Inset: examples of magnetization 
dynamics. 



where the rj are Gaussian ^-correlated noise: (ijq"\t) r ]l^ \t')) = ^a,a'^q,q'S(t — tf), and 
(r)o(t)r)o(t?)) = 2S(t — t'). We have now all the ingredients to compute P(m, r), the PDF of the 
instantaneous magnetization m at time r. It is given by, P(m,r) = (6(m — Tntot{T~))) {r)o,ri }• 
As previously [3,6] we introduce a reduced scaling variable for the magnetization. Here we are 
interested in fluctuations around the typical value of m, m = (m-tot) { nq } , where the average 
is performed over all noise except 770- Similarly we define the width of the distribution as 



a 2 = {TTLtot) {r]^} — m?! and the reduced magnetization /it c 



1X1)1 a. We stress that m 



we compare the 
with the results of 



and a are not the mean and variance of the instantaneous magnetization, as we have not 
averaged over the noise 770 ■ 

At this stage of our derivation, it is useful to numerically check the approximations and 
assumptions made so far, i.e. the perturbation expansion, the Langevin dynamics and the 
logic of the separation of time scales. In figure 1(a) we show data generated by integrating 
numerically equations 10} for a temperature below the mean-field critical temperature, at 
zero field. This shows that perturbation theory can indeed capture the first departure from 
Gaussian fluctuations, as the critical point is approached. In figure 
PDF obtained from Monte Carlo simulations at T*(L = 64) = T c x 0.94 l 
Langevin dynamics at T*(L = 65, A = 2.5) = T c mf x 0.95. As one can see the agreement is good 
and the data give an equally good fit to the BHP function. Differences appear in the wings 
of the distribution. The perturbation scheme overestimates the value of the PDF, illustrating 
the limit of its validity. However, it is clear that, at least in the case where the symmetry is 
effectively broken and where our criterion J5J is satisfied, the scheme captures the magnetic 
fluctuations of the Ising model to an excellent approximation. Our scheme has to be compare 
with mean-field treatment of Zheng [13] where a simple approximation or ansatz can also 
describe most of the features. Here however, the connection with the XY-model is revealed 
through a separation of the "fast" variables, described by spin waves like excitations, and 
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"slow" variables, representing the evolution of the global magnetization within the effective 
potential minima. 

Probability density function of the fluctuations. Using the integral representation of 
the Dirac function, P(m, t) can be written as a path integral over the noise r]o and rj q . As 
the equations on <j> q °^ are linear they can be integrated out if we assume that cf>o is slowly 
varying with time. This is true at low temperature: 4>o is the instantaneous magnetization 
within the mean-field approximation, and the fluctuations around this value are represented by 
the . The amplitude of the fluctuations of the global variable 4>o m equation are scaled 
by the factor 1/N coming from So, and therefore 4>o does not, at least at low temperature, 
venture far from the energy minima. In this case, the quantity G q , which depends on <po, can 
be considered as constant. That is, we can replace the propagator by its time averaged value 
during the period r. 

After some algebra, defining = 2 J2 qeSl [KqGq coth(r/G q )] fc /N k , one finally obtains 



P(m,r) = 
IIxY^oK^tot, r) 



V m (t) exp [~\J »?g(«)d«J Tl XY 



[</>0](Mtot,T), 



(5) 



+00 



dx / 

— exp ^ tot , 2 



9^ h \ \ n„ 



k=2 



The function IIxy is the PDF for the 2D-XY model in the spin wave approximation with a 
massive propagator, G~ x oc q 2 + 7\/(</> ) 2 [14]. Here IIxy depends on the temperature as the 
mass varies with temperature through <j>$. In the limit M — > 0, the temperature dependence 
disappears and IIxy becomes precisely the BHP function. 

In order to obtain the PDF for m to t, we would have to evaluate the last path integral over 
?7o . This integral is related to the non linear Langcvin equation (0J. At T = the dominant 
solution of the equation of motion is 4> °\t) = constant, which is a solution of 5So/64>o(u) = 0: 
the mode <J)q does not have any dynamics and it affects the PDF only by imposing a finite 
mass M(0q -'). Physically this means that the PDF for the Ising model at low temperature 
should be the same as that obtained for the 2D-XY model with a small magnetic field. This 



is exactly what we have observed in figure 1 (a) At finite temperature however non constant 
solutions exist: these are the instantons associated with the non- linear Langevin equation 
Expressing the path integral in © over T)o(t) as a path integral over <f>o(t), using equation J2J|, 
it follows that the time dependent solutions extremize the integral [4>o(u) + 5Sq/ '54>o(u)] 2 du. 
This leads to 



P(m,r) = / £><Mt)exp 



-\{S [Mr)]-S [MO)}) 



(6) 



exp 



<t> 2 (u) + 



SSn 



6<po(u) 



- 2 



du 



IIxy [<Ao](Mtot, 



The second order derivative of So, which comes from the Jacobian of the transformation, 
is of order TV compared with order N 2 for the other terms and is therefore negligible. To 
simplify the analysis, we assume periodic boundary conditions, </>o(0) = </>o( r )j and the first 
argument of the exponential in © vanishes. The non constant solutions <p k \u), k > 1 
verify </> (u) 2 /2 = Vo — V e s[4>o(u)], where Vb > is a constant of motion that depends on the 
instanton trajectory and 2V c tf[4>o(u)] = — [SSo/S(t)o(u)] 2 is an effective inverted potential. From 
(0 we define an effective action S c s whose (j)^ are the classical time dependent solutions: 
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Ses[(f>o(u)] = ^0o( w ) — K:ff [0o( u )]j dit. Expressing du as function of d</> , the constant 

Vo satisfies the equation J ± . = r. The sign is positive (resp. negative) when </>o is 

V 2(Vb — Ktf) 

increasing (resp. decreasing) with time. Generally, we assume that Vq tends exponentially to 
zero when r is large. For an action with two wells So = — a0o/2+0o/4, with a cx T c — T, we can 
evaluate this constant. Indeed, if we consider a trajectory going from <^o(0) = \fo> to (f>o = 
at some later fixed time, then going back to i/a at time t, this corresponds to solving the 

following equation: 2 J n , dl ^ Q = r. When r is large, Vo has to be small, and the main 

y/2(Vo — V ff) 

contributions from the previous integral come from the end points: V c ff ~ —2a 2 {4>o T V^) 2 
for </>o close to ±y/a, and V^g ~ — a 2 0g/2 for </>o close to 0. It is then easy to show that 
Vo oc exp(— ^/2ar) and it is hence legitimate to discard this constant for large time r. The 

(k) 

effective action S c g can be simplified in this case, and is equal to ASq , the sum of the energy 
barriers crossed by the instanton (fffi during the time r. Finally, the distribution 10 can be 
put into the following form 

P(m,r) ~ n X Y[^ 0) ](Mtot) + ^C fe cxp (-A5< fc) ) n XY [4 fe) ](Mtot,r), (7) 

k>l 

where the Cfe's arc coefficients related to the Gaussian fluctuations around the saddle point 
solutions (frb^ ■ These instantons restor the system symmetry. Below the critical region they 
appear only on exponentially large time scales, and can be neglected on the time scale of 
any observations. Near the critical point they appear more frequently and lead to a phase 
transition in the finite system. It is extremely challenging to exactly compute the contribution 
from the instantons [15] and to do so would, in any case, give only an approximate description 
of the true critical dynamics, so we do not attempt it here. Rather, the final expression is 
sufficient for our purposes as it exposes both the origin of generic critical fluctuations and of 
the dependence on universality classes through the generation of instanton like excitations. 

Interpretation and generalization. We have shown that at this level there are two 
distinct contributions to the PDF for magnetic fluctuations in equation @. The first is a 
Gaussian action coming directly from the perturbation expansion and as long as quadratic 
fluctuations arc present, such a term must appear. The limiting case of zero mass corresponds 
BHP distribution. The Ising model universality class does not appear in this term and in this 
sense it is superuniversal. The fluctuations described here are localized around the typical 
value of the magnetization and are sensitive to the local geometry around this point in phase 
space, not to the global structure of the phase space related to the universality class. The 
second term comes from the contribution of instantons that restore phase space symmetry 
and so is strongly dependent on the universality class. It is the analogue of the corrections 
computed exactly for the 2D-XY model within the spin wave approximation [9, 10] 
The existence of a point, or points, in the phase diagram where the PDF of the 2D Ising model 
is close to the BHP distribution [8] results from a compromise between these two terms: to be 
as close as possible to BHP, we have to reduce the mass in the Gaussian propagators of the first 
term. This is exactly what happens as one approaches T c and the instantaneous value of <f>o 
reduces. At the same time the contribution from instantons has to be small, which is not the 
case at T c : reducing <po increases the probability of inducing an instanton, taking the system 
from one local minimum to the other. This occurs just as the correlation length becomes of 
the order of the system size, hence the instanton contribution being small corresponds exactly 
to our criterion J2J being satisfied. The temperature T* studied in [8] is the point of best 
compromise. The application of a small magnetic field breaks the symmetry, increasing the 
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(k) 

barriers A5q . The contribution from instantons is then reduced and one can expect better 
agreement with BHP than in zero field. This is just what is observed. However, while the 
agreement with BHP can be excellent it is approximate and there will always be corrections 
to it. 

What other correlated systems show similar behaviour ? The previous analysis suggests 
that generic behaviour reminiscent of the D-dimensional Gaussian model should indeed be 
commonly observed in equilibrium correlated systems in D-dimensions. We show in fact that 
fluctuations of this generic type correspond to perturbative corrections to the central limit 
theorem and that in this sense such systems should be thought of as being weakly, rather 
than strongly correlated. It suggests also that one should expect variations. These variations 
can be large, coming from non-linear instanton-like objects, taking the distribution far from 



the generic asymmetric form shown in Fig 1(b) [8, 16], or they can be smaller, coming from 



the harmonic contribution itself. That is, the BHP function is not a miraculous unique 
function for all circumstances, rather distributions of the form of the first term in equation 
depend weakly on the boundary conditions [6, 17, 18], on the mass [14] and even on the 
temperature [9,10]. Of these parameters the strongest dependence is on dimension: we expect 
that 3D correlated equilibrium systems are related to the 3D Gaussian model, which has 
weakly asymmetric order parameter fluctuations and is not critical [6,19]. Interestingly, out 
of equilibrium three dimensional systems show results resembling the 2D-XY model [1,2,4], or 
the closely related Gumbel distribution [17] characteristic of a one-dimensional system with 
long range interactions not those of the 3D-XY model. These results seem to suggest the 
presence of a dimensional reduction in such systems that remains to be explained. 



It is a pleasure to thank S.T. Bramwell, S.T. Banks, F. Dclduc and P. Pujol for useful 
discussions and comments. 
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